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o : 

The Coulomb interaction effect and the surface effect are consistently taken into account 
£SJ ' in the hadron-quark mixed phase. These two finite-size effects greatly change the properties 

! i t of the mixed phase and restrict its density region. In particular, the charge screening effect 

^) . and the rearrangement of the charged particles are elucidated. Keeping the Gibbs conditions 

throughout the numerical procedure, we show the Maxwell construction effectively regain 
the physical meaning and the equation of state becomes similar to that given by the Maxwell 
04 , construction. 

<N . 

fNj ' §1. Introduction 

> : 

Deconfinement phase transition is believed to occur in hot and/or high-density 
£NJ ■ matter, while its mechanism has not been well understood yet. Many authors have 

studied this transition by model calculations or first-principle calculations like lattice 
\ QCD. 1 ) Nowadays it is widely accepted that quark matter exists in hot and/or 

high-density region like inner core of neutron stars. Static and dynamic properties 
of quark matter have been actively studied theoretically for quark-gluon plasma 
Q_i| (QGP), color superconductivity 2 )' 3 ) or magnetism. 4 )~ 6 ) Phenomenologically quark 

matter has been actively searched for in relativistic heavy-ion collisions (RHIC), 7 ) 
or in early universe and compact stars. 8 )' 9 ) 

Many theoretical calculations have suggested that the deconfinement phase tran- 
sition should be of first order in low temperature and high density area. 10 )' n ) There- 
fore we assume first-order phase transition in this study. We, hereafter, consider the 
phase transition from nuclear matter to three-flavor quark matter in neutron-star 
matter for simplicity. If the deconfinement transition is of first order, we may expect 
the mixed phase during the transition. The hadron-quark mixed phase has been 
considered during the hadronization in RHIC 12 )~ 14 ) or the boundary between quark 
matter and hadron matter in neutron stars. 15 - 1 

There is an issue about the mixed phase for the first-order phase transitions with 
more than one chemical potential. 16 ) We often use the Maxwell construction (MC) to 
derive the equation of state (EOS) in thermodynamic equilibrium, as in the water- 
vapor phase transition. In this case both phases consist of single particle species 
(H2O). However, if many particle species participate in the phase transition as in 
neutron-star matter, MC is no more an appropriate method. Before Glendenning 
first pointed out, 16 ) many people have applied MC to get EOS of the first-order 
phase transitions 17 ) -20 ) expected in neutron stars, such as pion or kaon condensation 
and the deconfinement transition. 

For the deconfinement transition in neutron-star matter, we consider quark de- 
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grees of freedom as well as hadrons and leptons. Accordingly we must introduce 
many chemical potentials for particle species, but the independent ones in this case 
are reduced to two, i.e. baryon-number chemical potential hb and charge chemical 
potential /j,q, due to beta-equilibrium and total charge neutrality. They are noth- 
ing but the neutron and electron chemical potentials, \i n and fi e , respectively. In 
the mixed-phase these chemical potentials should be spatially constant. When we 
naively apply MC to get EOS in thermodynamic equilibrium, we immediately notice 
that /xb is constant in the mixed-phase, while /i e is different in each phase because 
of the difference of the electron number in these phases. This is because MC uses 
EOS of bulk matter in each phase, which is of locally charge-neutral and uniform 
matter; many electrons are needed in hadron matter to cancel the positive charge 
of protons, while in quark matter total charge neutrality is almost fulfilled without 
electrons. Thus 

m§ = mI, (i-i) 

in MC, where superscripts "Q" and "H" denote the quark and hadron phase, re- 
spectively. Glendenning emphasized that we must use the Gibbs conditions (GC) in 
this case instead of MC, which relaxes the charge-neutrality condition to be glob- 
ally satisfied as a whole, not locally in each phase. 16 ** GC imposes the following 
conditions, 

,.Q _ ,.h Q _ H 

pQ — pH rpQ _ rpH (T2) 

He demonstrated a wide region of the mixed phase, where two phases have a net 
charge but totally charge-neutral: EOS thus obtained, different from that given by 
MC, never exhibits a constant-pressure region. He simply considered the mixed phase 
consisting of two bulk matters separated by a sharp boundary without any surface 
tension and the Coulomb interaction, which we call "bulk Gibbs" for convenience. 

"Bulk Gibbs" requires that each matter can have a net charge but total charge 
is neutral, 

fvpl + (1 - fv)pl = 0, (1-3) 

O H 

where fy means the volume fraction of quark matter in the mixed phase and "p^ " 
means charge density in each matter. Figure ^ shows the phase diagram in the hb 
- \i e plane. We can see that there is a discontinuous jump in fi e for the case of MC, 
while the curve given by "bulk Gibbs" smoothly connects uniform hadron matter 
and uniform quark matter; the mixed phase can appear in the wide /iB region in 
"bulk Gibbs", in contrast with MC. 21 ) 

However, this "bulk Gibbs" should be too simple to study the mixed phase, 
since we must consider non-uniform structures to be more realistic, instead of two 
bulk uniform matters; the mixed phase should have various geometrical structures, 
where both the number and charge densities are no more uniform. Then we have to 
take into account the finite-size effects like the surface and the Coulomb interaction 
energies. In this paper we study the structured mixed phase by treating the finite- 
size effects self-consistently. We show that the mixed phase should be narrow in 
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Fig. 1. Phase diagram in the /ib — Mo plane. There appears no region of the mixed phase by the 
calculation with MC, while a wide region of the mixed phase by the "bulk Gibbs" calculation. 

the //b space by the charge screening effect, and derive EOS for the deconfinement 
transition in neutron-star matter. We shall see it results in EOS being similar to the 
one given by MC. We also discuss the interplay of the Coulomb interaction effect 
and the surface effect in the context of the hadron-quark mixed phase. Preliminary 
results for the droplet case has been already reported in Ref. 28 - 1 

The plan of the paper is as follows. We briefly review the previous works in Sec. 
|2J Section |3] is devoted to our formalism and the numerical procedure. We show 
the results of our calculation and discuss the screening effect in Sec. Finally, a 
summary and concluding remarks are given in Sec. |SJ 



considered spherical quark droplets embedded in hadron matter by including the 
surface and the Coulomb energies. They introduced the surface tension and treated 
its strength as a free parameter because the surface tension at the hadron-quark in- 
terface has not been clearly understood. They pointed out that if the surface tension 
parameter a is large (a > 90 MeV/fm 2 ), the region of the mixed phase is largely 
limited or cannot exist. Subsequently Glendenning and Pei 15 ^ have suggested the 
"crystalline structures of the mixed phase" which have some geometrical structures, 
"droplet", "rod", "slab", "tube", and "bubble", assuming a small a. 15 )' 23 ) 

The finite-size effects are obvious in these calculations by observing energies. We 
may consider only a single cell, by dividing the whole space into equivalent Wigner- 
Seitz cells: the cell size is denoted by Rw and the size of the lump (droplet, rod, 
slab, tube or bubble) by R. Then the surface energy density is expressed in terms 
of the surface tension parameter a as 



§2. Brief review of the previous works 



Heiselberg et al. 22 ) studied a geometrical structure in the mixed phase: they 



es = 



fycrd 
R 



(2-1) 



where d denotes the dimensionality of each geometrical structure; d = 3 for droplet 
and bubble, d = 2 for rod and tube, and d = 1 for slab. The Coulomb energy density 
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reads 

(2-2) 

(d+2)~\ (2-3) 

where we simply assumed uniform density in each phase. When we minimize the 
sum of es and ec with respect to the size R for a given volume fraction fy, we can 
get the well known relation, 

es = 2ec- (2-4) 

This implies that an optimal size of the lump is determined by the balance of these 
finite-size effects. Eventually we can express the size-dependent energy density be- 
sides the bulk energy density: 22 ) 

Thus we can calculate the energy of any geometrically structured mixed phase with 
(|2-5jl by changing the parameter d. Many authors have taken this treatment for the 
mixed phase. 15 )' 21 )> 22 ) Note that the energy sum in Eq. (|2-5[) becomes larger as the 
surface tension gets stronger, while the relation Ea. (|2-4j) is always kept. 

However, this treatment is not a self-consistent, but a perturbative one, since the 
charge screening effect for the Coulomb potential or the rearrangement of charged- 
particles in the presence of the Coulomb interaction is completely discarded. We 
shall see that the Coulomb potential is never weak in the mixed phase, and thereby 
this treatment overestimates the Coulomb energy. The charge screening effect is 
included only if we introduce the Coulomb potential and consistently solve the Pois- 
son equation with other equations of motion for charged particles. Consequently it 
is a highly non-perturbative effect. Norsen and Reddy 24 ) have studied the Debye 
screening effect in the context of kaon condensation to see a large change of the 
charged-particle densities like kaons and protons. Maruyama et al. have numerically 
studied it in the context of liquid-gas phase transition at subnuclear densities, 25 ) 
where nuclear pastas can be regarded as geometrical structures in the mixed phase. 
Subsequently, they have also studied kaon condensation at high-densities, 26 ) where 
we have seen that kaonic pastas appear in the mixed phase. Through these works 
we have figured out the role of the Debye screening in the mixed phase. We have 
also studied the interplay of the Coulomb effect and the surface effect. 

Voskresensky et al. 27 ) explicitly studied the Debye screening effect for a few 
geometrical structures of the hadron-quark mixed phase. They have shown that 
the optimal value of the size of the structure cannot be obtained due to the charge 
screening even if the surface tension is not so strong. They called it as mechanical 
instability. It occurs because the Coulomb energy density is suppressed at larger 
size than the Debye screening length (cf. Eq. (|4-2|l ). They also suggested that the 
properties of the mixed phase become very similar to those given by MC, if the 



6 C = 2vre 2 (p c H h - p<i) R 2 Mfv] 



2(d-2) 



-i 



\dfy- 2/d ) +fv 
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charge screening effectively works. They also noted that the apparent violation of 
the Gibbs condition (Eq. can be remedied by including the Coulomb potential 

in a gauge-invariant way: the number of the charged particles is given by a gauge 
invariant combination of the chemical potential and the Coulomb potential, and 
thereby the number can be different in each phase for a constant charge chemical 
potential if the the Coulomb potential takes different values in both phases. However, 
they used a linear approximation to solve the Poisson equation analytically. 

If the Coulomb interaction effect is so important, it would be important to 
study it without recourse to any approximation. In this paper we numerically study 
the charge screening effect on the structured mixed phase during the deconfinement 
transition in neutron-star matter in a self-consistent way. Actually we shall see 
importance of non-linear effects included in the Poisson equation. 



§3. Self-consistent calculation 



nucleon phase 
quark phase 



3.1. Thermodynamic potential 

We consider the geometrically structured mixed phase (SMP) where one phase is 
embedded in the other phase with a certain geometrical form. 

We divide the whole space into equivalent charge- 
neutral Wigner-Seitz cells with a size Rw and a 
size of embedded phase R as illustrated in Fig. |2J 
Quark phase consists of u, d, s quarks and elec- 
tron. Hadron phase consists of proton, neutron 
and electron. We incorporate the MIT Bag model 
and assume the sharp boundary at the hadron- 
quark interface. We use density functional theory 
(DFT) and incorporate local density approxima- 
tion (LDA). 29 )< 3 °) 

We consider total thermodynamic potential 
(l? to tal) which consists of the hadron, quark and 
electron and the Coulomb interaction contribu- 
tions: 




nucleon quark 



Fig. 2. Wigner-Seitz approximation 
for the droplet case. R is the 
droplet radius and Rw the cell 
radius. 



^total — ^hadron ~i~ ^quark ~i~ 



(3-1) 



where we summarize the contributions of elec- 
trons and the Coulomb interaction as f2 cm because 
they are present in both phases. 

We briefly present the expressions of thermodynamic potentials. The details of 
the derivation of the expressions are given in Ref. 31 ) 

First, the Coulomb interaction energy is expressed in terms of particle densities, 



1 



iPi(r)QjPj{f) 



where i = u,d,s,p,n,e with Qi being the particle charge (Q 



(3-2) 

-e < for the 
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electron). Accordingly the Coulomb potential is denned as 



(3-3) 



where Vq is an arbitrary constant representing the gauge degree of freedom. We 
fix the gauge by a condition V(Rw) = in this paper (see Sec. 2.2). Operating a 
Laplacian V 2 on the Coulomb potential V(r), we automatically derive the Poisson 
equation. 

Therefore, the electron contribution and the Coulomb interaction energy (in 
both phases) are expressed as 



i 



8ire 2 
1 



8vre 2 



(W(r)) 2 + e e (Pe(r» - A*ePe(f) + V(r)p e (r) 



12vr 2 



(3-4) 



where e e (p e (r)) = ± — ^ 2 ' is the kinetic energy density of electron. 

Secondly, in the quark phase, u and d quarks are treated as massless particles 
and only s quark massive one, m s = 150 MeV. The kinetic energy of quark of flavor 
/ is simply expressed as 32 ) 



e /kin 



8vr 2 



mj [x f rjf (2xj + 1) - In (x f + Vf )] 



(3-5) 



where mj is the quark mass, Xf = ppf(r)/mf with Fermi momentum ppf(r) = 

(ir 2 pf(f))Ti and 77/ = yjl + xj. 

For the interaction energy, we take into account the leading order contribution 
coming from the one-gluon exchange interaction. Since the contribution from the 
Hartree term disappears due to the traceless property of the Gell-Mann matrix (A), 
the leading order contribution only comes from the Fock term, 

e/Fock = - ^mj I x) - I [x f T)f - In (x f + r] f )} 2 ^ . (3-6) 

Including this interaction, the quark contribution to the thermodynamic potential is 
expressed as 



^quark = + + + J dfB, 



(3-7) 



dr[e f ( Pf (f)) - p fPf {r) - N t V(r)p f (r)] , N t = % 



where, the energy density Cf(pf(r)) stands for ef^m + e/Fock of / quark, and B is 
the Bag constant. The Bag constant is taken as 120 MeV/fm 3 , and the QCD fine 
structure constant as a c = 0.4, which are also used by Heiselberg et al. 22 ) and in 
the previous work. 27 ) 
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Thirdly, we consider the hadron contribution. The thermodynamic potential for 
the non-relativistic nucleons becomes 

^hadron = E N - ^2 Pa J dfp a (f)-J dfV(f)p p (f), (3-8) 



a=p,n 



where En is the energy of the nucleons, 



En 



I 



df 



£ M 1 P2<f) + e P ot (p P (r-),p n (r)) 



a=p,n 



(3-9) 



Here we use the effective potential e po t(p p (r) , p n (f)) parametrized by the local den- 
sities for simplicity, 

(Pn(r) - p p (f)) 



e P ot(r) = s < 



+ (Pn(r) +p P (r))e bind 

2 



Po(0 

^ (Pn(r) + p p (r)) / p n (f) + pp N 
+ ^° 18 {—^ \ 

+ C sat (p n (r) + p p (r)) f Pn(0 + Pp(r1 - l) , (3-10) 
V Po J 

where So, Ko, eund, arid C sa t are adjustable parameters to reproduce the saturation 
properties of nuclear matter. 27 ) We consider beta equilibrium at the hadron-quark 
interface as well as in each phase: 

p u + p c = p d , 

Md = Ms; 

Pp + Pc = Pn = PB, 
Pn = Pu + 2/U d , 

p p = 2p u + p d . (3-11) 

The last relation can be derived from other four relations, so that there are left four 
independent conditions for chemical equilibrium. 

We get the equations of motion from **^ tal = ( fa = p u (f), pd(f), p s (f), 
p p (r), p n (f), p e (r),V(r) ): The Poisson equation then reads 



VV(r)=4^e 2 



2 1 1 

^Pu(r) - ~pd{r) - -p s {r) ) + P P (r) - pe{r) 



(3-12) 



Other equations of motion give nothing but the expressions of the chemical poten- 
tials, 

(3-14) 



as 



Pu 



i=u i d i s,e 

2a 
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Pd = ( 1+ ^ ) 7T3 P |(fO + -F(r) (3-15) 



. 2a c 

Ps = e Fs r + — 

37T 



PFs(^) - 3 In ' 



e Fs (r) V m s 



(3-16) 



with efs (0 = Y m s + Pfs 00 • On ^ ne other hand chemical potentials of nucleons and 
electrons are 

P 2 Fll , 2So(pn(r) -pp(r)) 

Pn — o 1 r e bind 

2m po 
+ Ko / pn(r)+p P (F) _ \ 2 + £o / Pn(f) + pp(f) _ 1 
6 V Po / 9 V Po 

+ 2C sat ^ + ^ (f) -C sat (3-17) 
Po 

„ _„ g|n(0, g|pg 45 (p B -2p p (r-)) 2 
Mp_/ " n 2m + 2m p V[r) 



p e ={M 2 Pe {r)Y +V{r). (3-18) 

We solve these equations of motion under GC. Important point is that the Coulomb 
potential V(r) is included in each expression in a proper way. The Coulomb potential 
is the function of charged-particle densities, and in turn densities are functions of 
the Coulomb potential. As a result, the Poisson equation becomes highly non-linear. 
Since it should be difficult to solve them analytically, we numerically solve them 
without any approximation. 

Once the geometrical structure is concerned, we have to take into account the 
surface tension as well at the interface of the hadron and quark phases. It may be 
connected with the confining mechanism and unfortunately we have no definite idea 
about how to incorporate it. Actually many authors have treated its strength as a 
free parameter and seen how the results are changed by its value. 15 )' 21 )' 22 ) Here we 
also follow this manner by introducing the surface tension parameter a to simulate 
the surface effect. One might be afraid that the surface tension will be modified once 
the Coulomb interaction is explicitly introduced. However, such modification might 
be rather small, as inferred from the previous result. 26 ) 

Note that we have to determine now eight variables, i.e., six chemical potentials, 
Pu> Pd, Psj Pp, Pnj Pe, and radii R and Rw- First, we fix R and Rw- Here we have 
four conditions due to (5 equilibrium ()3-ll|) . Therefore, once two chemical potentials 
/xb and /i e are given, we can determine other four chemical potentials, /x u , fid, Ps 
and fip. Next, we determine /i e by the global charge neutrality condition: 



fvpl + (1 - fv)Pch = 0, (3-19) 
geometrical structure. At this point fy is still fixed. 



where the volume fraction fy = (^jj^J , and d denotes the dimensionality of each 
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The pressure coming from the surface tension is given by 

P.-*§r, (3.20) 

where S is the area of the surface and Vq is the volume of the quark phase. Then we 
find the optimal value of R (Rw is fixed and thereby fy is changed by R) by using 
one of GC; 

pQ = p H + P (J . (3-21) 

The pressure in each phase pQ( H ) is given by the thermodynamic relation: pQ( H ) = 
— J7q( H )/Vq( H ), where ^q(h) is the thermodynamic potential in each phase and given 
by adding electron and the Coulomb interaction contributions to ^ qua k(hadron) i n Fqs. 
1)3-71) and (|3-8j) . Finally, we determine Rw by minimizing thermodynamic potential. 
Therefore once pb is given, all other values pn (i = u,d,s,p,e) and R, Rw can be 
obtained. 

Note that we keep GC throughout the numerical procedure. We will see later 
how the mixed phase would be changed by including finite-size effects keeping GC 
completely. Although MC is not rigidly correct as we have seen in Sec. 1, our results 
will show a similar behavior to those by MC as a result of including the finite-size 
effects. 

In numerical calculation, every point inside the cell is represented by a grid 
point (the number of grid points N gr ^ ~ 100). Equations of motion are solved by a 
relaxation method for a given baryon-number chemical potential under constraints 
of the global charge neutrality. 

3.2. Proper treatment of the Coulomb interaction 

With the Coulomb potential ()3-3)) and thermodynamic potentials ()3-l)) . the gauge 
invariance of our treatment can easily be seen as follows: varying the expression of 
chemical potentials 1)3-13)1 with respect to the Coulomb potential V(r), as is shown 
in the previous work, 27 ) we have 

dp ' 

Aij-rjy = Ni, AijBj k = Sik, (3-22) 
where matrices A and B are defined as 

dpidpj dpj 
From these equations, the gauge- invariance relation follows, 

We can understand that chemical potential is gauge variant from this relation. When 
the Coulomb potential is shifted by a constant value, V(f) ==>■ V{f) — Vq, the charge 
chemical potential should be also shifted as pi =>• pi + NiVq. To take into account 
the Coulomb interaction, we have to include V{r) in the gauge invariant way like in 
Eq. ()3-13)) . Note that the phase diagram in the pb — p e plane (see ,e.g., Fig. is not 
well-defined, since the charge chemical potential p e is not gauge invariant by itself. 
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§4. Numerical results 



We show the thermodynamic potential in Figs. |3] and 21 in uniform matter, 
hadron phase is thermodynamically favorable for lib < 1225 MeV and quark phase 
for /xb > 1225. Therefore we plot Su, difference of the thermodynamic potential 
density between the mixed phase and each uniform matter: 



5uj 



{^total 
Wtotal 



uniform 
H 

uniform 

Q 



/MB < 1225MeV, 
lib > 1225MeV, 



(4-1) 



where w tota i = Qotal/Vwi e * c - There we also depict two results for comparison: 
one is given by the "bulk Gibbs" calculation, where the finite-size effects are com- 
pletely discarded. The other is the thermodynamic potential given by a perturbative 
treatment of the Coulomb interaction, which is denoted by "no Coulomb" ; discard- 
ing the Coulomb potential V(r), we solve the equations of motion to get density 
profiles, then evaluate the Coulomb interaction energy by using the density profiles 
thus determined. We can see the screening effects by comparing this "no Coulomb" 
calculation with the self-consistent one denoted by "screening" . 8tu given by MC ap- 
pears as a point denoted by a circle in Figs. Eland where two conditions, = P H 
and /Xg = /Xg, are satisfied. On the other hand the mixed phase derived from "bulk 
Gibbs" appears in a wide region of zxb- Therefore, if the region of the mixed phase 
becomes narrower, it signals that the properties of the mixed phase become close to 
those of MC. One may clearly see that uxtotai becomes close to that given by MC due 
to the finite-size effects, the effects of the surface tension and the Coulomb interac- 
tion. The large increase of 5uj from the "bulk Gibbs" curve comes from the effects of 



> 

0J 



3 



10 



a= 40 MeV/fm 2 MC 




screening . 
"droplet; 
"rod 

slab 



tube 
bubble 



no Coulomhr 
■■ droplet! 
■■ rod 
■■ slab 
■■ tube 
■■bubble- 



1100 



1200 1300 
M-b [MeV] 



1400 



Fig. 3. Difference of the thermodynamic poten- 
tial density as a function of baryon-number 
chemical potential /xb for a = 40 MeV/fm 2 . 
If 8u> is negative, the mixed phase is a thermo- 
dynamically favorable state. MC determines 
one point of phase transition in uniform mat- 
ter, denoted as a circle in the [ib-Sui plane. 
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Fig. 4. Same as Fig. 3 for a = 60 MeV/fm 2 . 
The negative Su> region is narrower than the 
a = 40 MeV/fm 2 case. 



the surface tension and the Coulomb potential. Since the surface tension parameter 
is introduced by hand, we must carefully study the effects of the surface tension and 
the Coulomb interaction, separately. From the difference between the result given 
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by "no Coulomb" and that by "bulk Gibbs", we can roughly say that about 2/3 of 
the increase comes from the effect of the surface tension and 1 /3 from the Coulomb 
interaction (see Eq. (|2-4j) ). 

Comparing the result of self-consistent calculation with that of "no Coulomb" , 
we can see that the change of energy caused by the screening effect is not so large, but 
still the same order of magnitude as that given by the surface effect. If the surface 
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p d xl00 
p s xl00 



> 100 
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Fig. 5. Density profiles in the droplet phase Fig. 6. Density profiles and the Coulomb poten- 
given by "no Coulomb" for = 1189 MeV tial given by the self-consistent calculation for 
and a = 60 MeV/fm 2 . They are uniform in the same parameter set as Fig. 5. R = 7.7 fm 
each phase. R = 7.2 fm and Rw = 12.8 fm. and Rw = 18.9 fm. 



tension is stronger, the relative importance of the screening effect becomes smaller 
and the effect of the surface tension becomes more dominant, as is seen in Fig. 4. 
To be more realistic we have to take 
into account the modification of the 1 
surface tension as the structure size 
changes. Though we cannot clearly 
say how the surface tension is af- 
fected, we may infer from the previous 'g 
study that it is not so large. Although '— ' 
the charge screening has not so large 
effects on bulk properties of the mat- 
ter, we shall see that it is remarkable — \ 
for the charged particles to change the 
properties of the mixed phase. 

The screening effect induces the Fig. 
rearrangement of the charged parti- 
cles. We can see this screening ef- 
fect by comparing Fig. \5\ with Fig. H3 
The quark phase is negatively charged 
and the hadron phase is positively 
charged. The negatively charged par- 
ticles in the quark phase such as d, s, 
e and the positively charged particle 





p h screening 
p„ screening 

p h no Coulomb 

p q no Coulomb 








10 
r [fm] 



20 



7. Local-charge densities for the "no 
Coulomb" case and the case of the self- 
consistent with the screening effect. For 
"no Coulomb", charge density of each phase 
is constant over the region. The absolute 
value of the charge density is larger than that 
given by the self-consistent calculation in each 
phase. In the hadron phase, the charge den- 
sity becomes almost vanished near the cell 
boundary r = Rw- 
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in the hadron phase p are attracted 

toward the boundary. On the contrary the positively charged particle in the quark 
phase u and negatively charged particle in the hadron phase e are repelled from the 
boundary. The charge screening effect also reduces the net charge in each phase. In 
Fig. [7[ we show the local charge densities of the two cases shown in Figs. and El 
The change of the number of charged particles due to the screening is as follows: In 
the quark phase, the numbers of d and s quarks and electrons decrease, while the 
number of u quark increases. In the hadron phase, on the other hand, the proton 
number should decrease and the electron number should increase. Consequently the 
local charge decreases in the both phases. In Fig. [7|we can see that the core region 
of the droplet tends to be charge-neutral and near the boundary of the Wigner-Seitz 
cell is almost charge-neutral. 
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Fig. 8. Lump and cell radii given by the "no 
Coulomb" calculation. The Debye screening 
length is also depicted for comparison. R is 
thick-solid line and Rw is thick-dashed line. 
We can see the size of the structure becomes 
less than the Debye screening length. 
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Fig. 9. Same as Fig. 8 given by the self- 
consistent calculation with the screening ef- 
fect. The size of the structure becomes larger 
than that given by "no Coulomb" , and conse- 
quently exceeds the Debye screening length. 



In Figs. IBlandElwe present 
the lump and cell radii for each 
density. As we have shown 
in the previous paper, 27 ) the 
Coulomb energy is suppressed 
for larger R by the screening ef- 
fect. The R dependence of the 
total thermodynamic potential 
comes from the contributions 
of the surface tension and the 
Coulomb interaction: the opti- 
mal radius giving the minimum 
of the thermodynamic poten- 
tial is then determined by the 
balance between two contribu- 



Screening length 





Fig. 10. Schematic graphs of the droplet size and the 
Debye screening length. Right figure shows the case 
of the self-consistent calculation with the screening 
effect and left figure "no Coulomb". 
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tions, since the former gives a decreasing function, while the latter an increasing 
one. If the Coulomb energy is suppressed, the minimum of the thermodynamic po- 
tential is shifted to larger radius. As a result the size of the embedded phase (R) and 
the cell size (Rw) become large. In Ref. 27 ) they demonstrated that the minimum 
disappears for a large value of the surface tension parameter: the structure becomes 
mechanically unstable in this case. We cannot show it directly in our framework 
because such unstable solutions are automatically excluded during the numerical 
procedure, while we can see its tendency in Figs. |H1 and R and Rw get larger by 
the screening effect. 

We also see the relation between the size of the geometrical structure and the 
Debye screening length. The Debye screening length appears in the linearized Poisson 
equation and is then given as 

Kr-x>(^> <*r-*(3f). ^-^m 

(4-2) 

where (pf 1 ) stands for the averaged density in quark phase, (p c ^) is proton number 
averaged density in the hadron phase and (pf 1 ) is the electron charge density averaged 
inside the cell. It gives a rough measure for the screening effect: At a distance larger 
than the Debye screening length, the Coulomb interaction is effectively suppressed. 

In Fig. El we show sizes of geometrical structure for "no Coulomb" case. If we 
ignore the screening effect, the size of the embedded phase is comparable or smaller 
than the corresponding quark Debye screening length \ q D (Fig. I1U|). This may mean 
that the Debye screening is not so important. Actually, many authors have neglected 
the screening effect due to this argument. 21 )' 22 ) In Fig. 01 however, we see that the 
size of the embedded phase can be larger than X q D (Fig. ITO]) in the self-consistent 
calculation. We can also see the similar situation about Rw and X e D . This means 
that the screening has important effects in this mixed phase. We cannot expect such 
a effect without solving the Poisson equation because of the non-linearity. We show 
the EOS in Figs. ^2 and 1121 The pressure of the mixed phase becomes similar to 
that given by MC due to the screening effect. 

We have used a fixed surface tension parameter in the present study. Surface 
tension is a very difficult problem because it should be self-consistent with the two 
phases of matter, quark and hadron. Lattice QCD, based on the first principle, 
would be the most reliable theory. It predicts that the surface tension can be 10-100 
MeV/fm 2 . 33 )' 34 ) Although this range is for high temperature, our choice is within 
it. Moreover, other model calculations of the surface tension 35 ) -38 ) are similar to our 
choice. Although we cannot conclude that MC is perfectly correct, we can say that 
the results obtained by the "no Coulomb" calculation, which many authors have 
used, have to be checked again taking into account the finite-size effects. 

Let us consider some implication of these results for neutron star phenomena. 
Glendenning 15 ) suggested many SMP appear in the core region by using "bulk 
Gibbs": the mixed phase should appear for several kilometers. However we can 
say that the region of SMP should be narrow in the //b space and EOS is more sim- 
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Fig. 11. Pressure as a function of baryon- 
number density given by the "no Coulomb" 
calculation for a = 40MeV/fm 2 . The results 
given by "bulk Gibbs" and MC are also pre- 
sented for comparison. 



Fig. 12. Same as Fig. 1111 given by the self- 
consistent calculation with the screening ef- 
fect. 



ilar to that of MC due to the finite-size effects. These results correspond to recent 
other calculations. Bejger et al. 39 ) have examined the relation between the mixed 
phase and glitch phenomena, and shown that the mixed phase should be narrow if 
the glitch is generated by the mixed phase in the inner core. On the other hand the 
gravitational wave asks for density discontinuity in the core region. 40 ) These studies 
support our result. 

§5. Summary and concluding remarks 

We have numerically studied the charge screening effect in the hadron-quark 
mixed phase, by fully including the non-linear effects in the Poisson equation. Com- 
paring the results with those given by "no Coulomb" calculation, we have elucidated 
the screening effect. The density profiles of the charged particles are much mod- 
ified by the screening effect, while the thermodynamic potential is not so much 
affected; the charge rearrangement induced by the screening effect tends to make 
the net charge smaller in each phase. Consequently the system tends to be locally 
charge-neutral, which suggests that MC is effectively justified even if it is thermody- 
namically incorrect. In this context, it would be interesting to refer to the work by 
Heiselberg, 41 ) who studied the screening effect on a quark droplet (strangelet) in the 
vacuum, and suggested the importance of the rearrangement of charged particles. 

We have seen that thermodynamic quantities such as thermodynamic potential 
and pressure become close to those derived from MC by the screening effect, which 
also suggests that MC is effectively justified due to the screening effect. As another 
case of more than one chemical potential system, kaon condensation has been also 
studied 25 - 1 and the results are similar to those in the present study. Thus the impor- 
tance of the screening effect should be a common feature for the first-order phase 
transitions in high-density matter. 

We have included the surface tension at the hadron-quark interface, while its 
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definite value is not clear at present. There are also many estimations for the surface 
tension at the hadron-quark interface in lattice QCD, 33 )' 34 ) in shell-model calcula- 
tions 35 ^ 37 ^ and in model calculations based on the dual-Ginzburg Landau theory. 38 ) 
Our parameter is in that reasonable range. 

We have considered some implications of our results for neutron star phenomena. 
The screening effect would restrict the allowed SMP region in neutron stars, in 
contrast with a wide region given by "bulk Gibbs". 15 )' 23 ) It could be said that they 
should change the bulk property of neutron stars, especially the structure of the core 
region. 

Compact stars have the strong magnetic field and its origin is not well under- 
stood. One possibility is that it comes from quark matter in the core. 4 ) -6 ) Therefore 
it should be interesting to include the magnetic field contribution in our formalism. 

We have assumed zero temperature here. It would be much interesting to include 
the finite-temperature effect. Then it is possible to draw the phase diagram in the //b 
- T plane and we can study the properties of the deconfinement phase transition; our 
study may be extended to treat the mixed phase to appear during the hadronization 
of QGP in the nucleus-nucleus collisions and supernova explosions. 

In this study we have used a simple model for quark matter to figure out the 
finite-size effects in the SMP. However, it has been suggested that the color super- 
conductivity would be a ground state of quark matter. 2 )> 21 ) Hence we will include it 
in a further study. The hadron phase should be also treated more realistically; for 
example, we should include the hyperons or kaons in hadron matter. In the recent 
studies the mixed phase has been also studied 42 )' 43 ) in the context of various phases 
in the color superconducting phase. 
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